Drug-target identification in COVID-19 disease mechanisms using computational systems biology approaches.

Introduction: The COVID-19 Disease Map project is a large-scale community effort uniting 277 scientists from 130 Institutions around the globe. We use high-quality, mechanistic content describing SARS-CoV-2-host interactions and develop interoperable bioinformatic pipelines for novel target identification and drug repurposing. Methods: Extensive community work allowed an impressive step forward in building interfaces between Systems Biology tools and platforms. Our framework can link biomolecules from omics data analysis and computational modelling to dysregulated pathways in a cell-, tissue- or patient-specific manner. Drug repurposing using text mining and AI-assisted analysis identified potential drugs, chemicals and microRNAs that could target the identified key factors. Results: Results revealed drugs already tested for anti-COVID-19 efficacy, providing a mechanistic context for their mode of action, and drugs already in clinical trials for treating other diseases, never tested against COVID-19. Discussion: The key advance is that the proposed framework is versatile and expandable, offering a significant upgrade in the arsenal for virus-host interactions and other complex pathologies.

COVID19 Disease Map, a computational knowledge repository of virus-host interaction mechanisms.

We need to effectively combine the knowledge from surging literature with complex datasets to propose mechanistic models of SARS-CoV-2 infection, improving data interpretation and predicting key targets of intervention. Here, we describe a large-scale community effort to build an open access, interoperable and computable repository of COVID-19 molecular mechanisms. The COVID-19 Disease Map (C19DMap) is a graphical, interactive representation of disease-relevant molecular mechanisms linking many knowledge sources. Notably, it is a computational resource for graph-based analyses and disease modelling. To this end, we established a framework of tools, platforms and guidelines necessary for a multifaceted community of biocurators, domain experts, bioinformaticians and computational biologists. The diagrams of the C19DMap, curated from the literature, are integrated with relevant interaction and text mining databases. We demonstrate the application of network analysis and modelling approaches by concrete examples to highlight new testable hypotheses. This framework helps to find signatures of SARS-CoV-2 predisposition, treatment response or prioritisation of drug candidates. Such an approach may help deal with new waves of COVID-19 or similar pandemics in the long-term perspective.

Gene- and pathway-level analyses of iCOGS variants highlight novel signaling pathways underlying familial breast cancer susceptibility.

Single-nucleotide polymorphisms (SNPs) in over 180 loci have been associated with breast cancer (BC) through genome-wide association studies involving mostly unselected population-based case-control series. Some of them modify BC risk of women carrying a BRCA1 or BRCA2 (BRCA1/2) mutation and may also explain BC risk variability in BC-prone families with no BRCA1/2 mutation. Here, we assessed the contribution of SNPs of the iCOGS array in GENESIS consisting of BC cases with no BRCA1/2 mutation and a sister with BC, and population controls. Genotyping data were available for 1281 index cases, 731 sisters with BC, 457 unaffected sisters and 1272 controls. In addition to the standard SNP-level analysis using index cases and controls, we performed pedigree-based association tests to capture transmission information in the sibships. We also performed gene- and pathway-level analyses to maximize the power to detect associations with lower-frequency SNPs or those with modest effect sizes. While SNP-level analyses identified 18 loci, gene-level analyses identified 112 genes. Furthermore, 31 Kyoto Encyclopedia of Genes and Genomes and 7 Atlas of Cancer Signaling Network pathways were highlighted (false discovery rate of 5%). Using results from the "index case-control" analysis, we built pathway-derived polygenic risk scores (PRS) and assessed their performance in the population-based CECILE study and in a data set composed of GENESIS-affected sisters and CECILE controls. Although these PRS had poor predictive value in the general population, they performed better than a PRS built using our SNP-level findings, and we found that the joint effect of family history and PRS needs to be considered in risk prediction models.

Author Correction: COVID-19 Disease Map, building a computational repository of SARS-CoV-2 virus-host interaction mechanisms.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Comprehensive Map of the Regulated Cell Death Signaling Network: A Powerful Analytical Tool for Studying Diseases.

The processes leading to, or avoiding cell death are widely studied, because of their frequent perturbation in various diseases. Cell death occurs in three highly interconnected steps: Initiation, signaling and execution. We used a systems biology approach to gather information about all known modes of regulated cell death (RCD). Based on the experimental data retrieved from literature by manual curation, we graphically depicted the biological processes involved in RCD in the form of a seamless comprehensive signaling network map. The molecular mechanisms of each RCD mode are represented in detail. The RCD network map is divided into 26 functional modules that can be visualized contextually in the whole seamless network, as well as in individual diagrams. The resource is freely available and accessible via several web platforms for map navigation, data integration, and analysis. The RCD network map was employed for interpreting the functional differences in cell death regulation between Alzheimer's disease and non-small cell lung cancer based on gene expression data that allowed emphasizing the molecular mechanisms underlying the inverse comorbidity between the two pathologies. In addition, the map was used for the analysis of genomic and transcriptomic data from ovarian cancer patients that provided RCD map-based signatures of four distinct tumor subtypes and highlighted the difference in regulations of cell death molecular mechanisms.

A multiscale signalling network map of innate immune response in cancer reveals cell heterogeneity signatures.

The lack of integrated resources depicting the complexity of the innate immune response in cancer represents a bottleneck for high-throughput data interpretation. To address this challenge, we perform a systematic manual literature mining of molecular mechanisms governing the innate immune response in cancer and represent it as a signalling network map. The cell-type specific signalling maps of macrophages, dendritic cells, myeloid-derived suppressor cells and natural killers are constructed and integrated into a comprehensive meta map of the innate immune response in cancer. The meta-map contains 1466 chemical species as nodes connected by 1084 biochemical reactions, and it is supported by information from 820 articles. The resource helps to interpret single cell RNA-Seq data from macrophages and natural killer cells in metastatic melanoma that reveal different anti- or pro-tumor sub-populations within each cell type. Here, we report a new open source analytic platform that supports data visualisation and interpretation of tumour microenvironment activity in cancer.

The Phosphatase PRL-3 Is Involved in Key Steps of Cancer Metastasis.

PRL-3 belongs to the PRL phosphatase family. Its physiological role remains unclear, but many studies have identified that PRL-3 is a marker of cancer progression and shown it to be associated with metastasis. Evidence implicating PRL-3 in various elements of the metastatic process, such as the cell cycle, survival, angiogenesis, adhesion, cytoskeleton remodeling, EMT, motility and invasion, has been reported. Furthermore, several molecules acting as direct or indirect substrates have been identified. However, this information was obtained in many different studies, and it remains difficult to see the larger picture. We therefore systematically collected the published information together and used it to develop a comprehensive signaling network map. By analyzing this network map, we were able to retrieve the signaling pathways via which PRL-3 governs the key steps of the metastatic process in cancer. In this review, we summarize current knowledge of the role of PRL-3 in cancer and the molecular mechanisms involved. We also provide the web-based open-source PRL-3 signaling network map, for use in further studies.

Metabolic and signalling network maps integration: application to cross-talk studies and omics data analysis in cancer.

BACKGROUND: The interplay between metabolic processes and signalling pathways remains poorly understood. Global, detailed and comprehensive reconstructions of human metabolism and signalling pathways exist in the form of molecular maps, but they have never been integrated together. We aim at filling in this gap by integrating of both signalling and metabolic pathways allowing a visual exploration of multi-level omics data and study of cross-regulatory circuits between these processes in health and in disease. RESULTS: We combined two comprehensive manually curated network maps. Atlas of Cancer Signalling Network (ACSN), containing mechanisms frequently implicated in cancer; and ReconMap 2.0, a comprehensive reconstruction of human metabolic network. We linked ACSN and ReconMap 2.0 maps via common players and represented the two maps as interconnected layers using the NaviCell platform for maps exploration ( https://navicell.curie.fr/pages/maps_ReconMap%202.html ). In addition, proteins catalysing metabolic reactions in ReconMap 2.0 were not previously visually represented on the map canvas. This precluded visualisation of omics data in the context of ReconMap 2.0. We suggested a solution for displaying protein nodes on the ReconMap 2.0 map in the vicinity of the corresponding reaction or process nodes. This permits multi-omics data visualisation in the context of both map layers. Exploration and shuttling between the two map layers is possible using Google Maps-like features of NaviCell. The integrated networks ACSN-ReconMap 2.0 are accessible online and allows data visualisation through various modes such as markers, heat maps, bar-plots, glyphs and map staining. The integrated networks were applied for comparison of immunoreactive and proliferative ovarian cancer subtypes using transcriptomic, copy number and mutation multi-omics data. A certain number of metabolic and signalling processes specifically deregulated in each of the ovarian cancer sub-types were identified. CONCLUSIONS: As knowledge evolves and new omics data becomes more heterogeneous, gathering together existing domains of biology under common platforms is essential. We believe that an integrated ACSN-ReconMap 2.0 networks will help in understanding various disease mechanisms and discovery of new interactions at the intersection of cell signalling and metabolism. In addition, the successful integration of metabolic and signalling networks allows broader systems biology approach application for data interpretation and retrieval of intervention points to tackle simultaneously the key players coordinating signalling and metabolism in human diseases.

Creation and analysis of biochemical constraint-based models using the COBRA Toolbox v.3.0.

Constraint-based reconstruction and analysis (COBRA) provides a molecular mechanistic framework for integrative analysis of experimental molecular systems biology data and quantitative prediction of physicochemically and biochemically feasible phenotypic states. The COBRA Toolbox is a comprehensive desktop software suite of interoperable COBRA methods. It has found widespread application in biology, biomedicine, and biotechnology because its functions can be flexibly combined to implement tailored COBRA protocols for any biochemical network. This protocol is an update to the COBRA Toolbox v.1.0 and v.2.0. Version 3.0 includes new methods for quality-controlled reconstruction, modeling, topological analysis, strain and experimental design, and network visualization, as well as network integration of chemoinformatic, metabolomic, transcriptomic, proteomic, and thermochemical data. New multi-lingual code integration also enables an expansion in COBRA application scope via high-precision, high-performance, and nonlinear numerical optimization solvers for multi-scale, multi-cellular, and reaction kinetic modeling, respectively. This protocol provides an overview of all these new features and can be adapted to generate and analyze constraint-based models in a wide variety of scenarios. The COBRA Toolbox v.3.0 provides an unparalleled depth of COBRA methods.

Application of Atlas of Cancer Signalling Network in preclinical studies.

Cancer initiation and progression are associated with multiple molecular mechanisms. The knowledge of these mechanisms is expanding and should be converted into guidelines for tackling the disease. Here, we discuss the formalization of biological knowledge into a comprehensive resource: the Atlas of Cancer Signalling Network (ACSN) and the Google Maps-based tool NaviCell, which supports map navigation. The application of ACSN for omics data visualization, in the context of signalling maps, is possible via the NaviCell Web Service module and through the NaviCom tool. It allows generation of network-based molecular portraits of cancer using multilevel omics data. We review how these resources and tools are applied for cancer preclinical studies. Structural analysis of the maps together with omics data helps to rationalize the synergistic effects of drugs and allows design of complex disease stage-specific druggable interventions. The use of ACSN modules and maps as signatures of biological functions can help in cancer data analysis and interpretation. In addition, they empowered finding of associations between perturbations in particular molecular mechanisms and the risk to develop a specific type of cancer. These approaches are helpful, among others, to study the interplay between molecular mechanisms of cancer. It opens an opportunity to decipher how gene interactions govern the hallmarks of cancer in specific contexts. We discuss a perspective to develop a flexible methodology and a pipeline to enable systematic omics data analysis in the context of signalling network maps, for stratifying patients and suggesting interventions points and drug repositioning in cancer and other diseases.

Community-driven roadmap for integrated disease maps.

The Disease Maps Project builds on a network of scientific and clinical groups that exchange best practices, share information and develop systems biomedicine tools. The project aims for an integrated, highly curated and user-friendly platform for disease-related knowledge. The primary focus of disease maps is on interconnected signaling, metabolic and gene regulatory network pathways represented in standard formats. The involvement of domain experts ensures that the key disease hallmarks are covered and relevant, up-to-date knowledge is adequately represented. Expert-curated and computer readable, disease maps may serve as a compendium of knowledge, allow for data-supported hypothesis generation or serve as a scaffold for the generation of predictive mathematical models. This article summarizes the 2nd Disease Maps Community meeting, highlighting its important topics and outcomes. We outline milestones on the roadmap for the future development of disease maps, including creating and maintaining standardized disease maps; sharing parts of maps that encode common human disease mechanisms; providing technical solutions for complexity management of maps; and Web tools for in-depth exploration of such maps. A dedicated discussion was focused on mathematical modeling approaches, as one of the main goals of disease map development is the generation of mathematically interpretable representations to predict disease comorbidity or drug response and to suggest drug repositioning, altogether supporting clinical decisions.

The Virtual Metabolic Human database: integrating human and gut microbiome metabolism with nutrition and disease.

A multitude of factors contribute to complex diseases and can be measured with 'omics' methods. Databases facilitate data interpretation for underlying mechanisms. Here, we describe the Virtual Metabolic Human (VMH, www.vmh.life) database encapsulating current knowledge of human metabolism within five interlinked resources 'Human metabolism', 'Gut microbiome', 'Disease', 'Nutrition', and 'ReconMaps'. The VMH captures 5180 unique metabolites, 17 730 unique reactions, 3695 human genes, 255 Mendelian diseases, 818 microbes, 632 685 microbial genes and 8790 food items. The VMH's unique features are (i) the hosting of the metabolic reconstructions of human and gut microbes amenable for metabolic modeling; (ii) seven human metabolic maps for data visualization; (iii) a nutrition designer; (iv) a user-friendly webpage and application-programming interface to access its content; (v) user feedback option for community engagement and (vi) the connection of its entities to 57 other web resources. The VMH represents a novel, interdisciplinary database for data interpretation and hypothesis generation to the biomedical community.

Systems medicine disease maps: community-driven comprehensive representation of disease mechanisms.

The development of computational approaches in systems biology has reached a state of maturity that allows their transition to systems medicine. Despite this progress, intuitive visualisation and context-dependent knowledge representation still present a major bottleneck. In this paper, we describe the Disease Maps Project, an effort towards a community-driven computationally readable comprehensive representation of disease mechanisms. We outline the key principles and the framework required for the success of this initiative, including use of best practices, standards and protocols. We apply a modular approach to ensure efficient sharing and reuse of resources for projects dedicated to specific diseases. Community-wide use of disease maps will accelerate the conduct of biomedical research and lead to new disease ontologies defined from mechanism-based disease endotypes rather than phenotypes.

Signalling maps in cancer research: construction and data analysis.

Generation and usage of high-quality molecular signalling network maps can be augmented by standardizing notations, establishing curation workflows and application of computational biology methods to exploit the knowledge contained in the maps. In this manuscript, we summarize the major aims and challenges of assembling information in the form of comprehensive maps of molecular interactions. Mainly, we share our experience gained while creating the Atlas of Cancer Signalling Network. In the step-by-step procedure, we describe the map construction process and suggest solutions for map complexity management by introducing a hierarchical modular map structure. In addition, we describe the NaviCell platform, a computational technology using Google Maps API to explore comprehensive molecular maps similar to geographical maps and explain the advantages of semantic zooming principles for map navigation. We also provide the outline to prepare signalling network maps for navigation using the NaviCell platform. Finally, several examples of cancer high-throughput data analysis and visualization in the context of comprehensive signalling maps are presented.

Fibroblast Heterogeneity and Immunosuppressive Environment in Human Breast Cancer.

Carcinoma-associated fibroblasts (CAF) are key players in the tumor microenvironment. Here, we characterize four CAF subsets in breast cancer with distinct properties and levels of activation. Two myofibroblastic subsets (CAF-S1, CAF-S4) accumulate differentially in triple-negative breast cancers (TNBC). CAF-S1 fibroblasts promote an immunosuppressive environment through a multi-step mechanism. By secreting CXCL12, CAF-S1 attracts CD4+CD25+ T lymphocytes and retains them by OX40L, PD-L2, and JAM2. Moreover, CAF-S1 increases T lymphocyte survival and promotes their differentiation into CD25HighFOXP3High, through B7H3, CD73, and DPP4. Finally, in contrast to CAF-S4, CAF-S1 enhances the regulatory T cell capacity to inhibit T effector proliferation. These data are consistent with FOXP3+ T lymphocyte accumulation in CAF-S1-enriched TNBC and show how a CAF subset contributes to immunosuppression.

NaviCom: a web application to create interactive molecular network portraits using multi-level omics data.

Human diseases such as cancer are routinely characterized by high-throughput molecular technologies, and multi-level omics data are accumulated in public databases at increasing rate. Retrieval and visualization of these data in the context of molecular network maps can provide insights into the pattern of regulation of molecular functions reflected by an omics profile. In order to make this task easy, we developed NaviCom, a Python package and web platform for visualization of multi-level omics data on top of biological network maps. NaviCom is bridging the gap between cBioPortal, the most used resource of large-scale cancer omics data and NaviCell, a data visualization web service that contains several molecular network map collections. NaviCom proposes several standardized modes of data display on top of molecular network maps, allowing addressing specific biological questions. We illustrate how users can easily create interactive network-based cancer molecular portraits via NaviCom web interface using the maps of Atlas of Cancer Signalling Network (ACSN) and other maps. Analysis of these molecular portraits can help in formulating a scientific hypothesis on the molecular mechanisms deregulated in the studied disease. Database URL: NaviCom is available at https://navicom.curie.fr.

Drug-Driven Synthetic Lethality: Bypassing Tumor Cell Genetics with a Combination of AsiDNA and PARP Inhibitors.

Purpose: Cancer treatments using tumor defects in DNA repair pathways have shown promising results but are restricted to small subpopulations of patients. The most advanced drugs in this field are PARP inhibitors (PARPi), which trigger synthetic lethality in tumors with homologous recombination (HR) deficiency. Using AsiDNA, an inhibitor of HR and nonhomologous end joining, together with PARPi should allow bypassing the genetic restriction for PARPi efficacy.Experimental Design: We characterized the DNA repair inhibition activity of PARPi (olaparib) and AsiDNA by monitoring repair foci formation and DNA damage. We analyzed the cell survival to standalone and combined treatments of 21 tumor cells and three nontumor cells. In 12 breast cancer (BC) cell lines, correlation with sensitivity to each drug and transcriptome were statistically analyzed to identify resistance pathways.Results: Molecular analyses demonstrate that olaparib and AsiDNA respectively prevent recruitment of XRCC1 and RAD51/53BP1 repair enzymes to damage sites. Combination of both drugs increases the accumulation of unrepaired damage resulting in an increase of cell death in all tumor cells. In contrast, nontumor cells do not show an increase of DNA damage nor lethality. Analysis of multilevel omics data from BC cells highlighted different DNA repair and cell-cycle molecular profiles associated with resistance to AsiDNA or olaparib, rationalizing combined treatment. Treatment synergy was also confirmed with six other PARPi in development.Conclusions: Our results highlight the therapeutic interest of combining AsiDNA and PARPi to recapitulate synthetic lethality in all tumors independently of their HR status. Clin Cancer Res; 23(4); 1001-11. ©2016 AACR.

Network-based approaches for drug response prediction and targeted therapy development in cancer.

Signaling pathways implicated in cancer create a complex network with numerous regulatory loops and redundant pathways. This complexity explains frequent failure of one-drug-one-target paradigm of treatment, resulting in drug resistance in patients. To overcome the robustness of cell signaling network, cancer treatment should be extended to a combination therapy approach. Integrating and analyzing patient high-throughput data together with the information about biological signaling machinery may help deciphering molecular patterns specific to each patient and finding the best combinations of candidates for therapeutic targeting. We review state of the art in the field of targeted cancer medicine from the computational systems biology perspective. We summarize major signaling network resources and describe their characteristics with respect to applicability for drug response prediction and intervention targets suggestion. Thus discuss methods for prediction of drug sensitivity and intervention combinations using signaling networks together with high-throughput data. Gradual integration of these approaches into clinical routine will improve prediction of response to standard treatments and adjustment of intervention schemes.